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ABSTRACT Results of a-particle observations performed by the European satellite HEOS 1 , in the period 
from December 9, 1968, to April 13, 1969, and from September 6, 1969, to April 15, 
1970, are presented. The average bulk velocities of protons Vp and a-p articles Re appear 
to be equal; however, due to an instrumental bias, the possibility of V a being lower than 
Vp cannot be ruled out. Comparison with observations of Vela 3 and Explorer 34 
satellites gives evidence of a dependence of helium abundance on the solar cycle. The 
problem of the stability of differences between the bulk velocities of protons and 
at-p articles is investigated. The behavior of a-particles through interplanetary shock waves 
is illustrated in connection with magnetic field measurements. 


INTRODUCTION 

Initial results on the observations of the helium compo- 
nent of the solar wind, performed by the ESRO satellite 
HEOS 1, are given by Formisano et al [1970a] and 
Formisano et al [19706] . Here we summarize the 
results obtained during the period from December 9, 
1968, to April 13, 1969, and from September 6, 1969, 
to April 15, 1970. Observations of three interplanetary 
shock waves are presented, using HEOS 1 plasma and 
magnetic field data. Details of the satellite orbit, the 
particle detector, and the measurement routine are given 
by Bonetti et.al [1969] , 

HEOS 1 was launched on December 5, 1968, on a 
highly eccentric orbit. The satellite is spin stabilized; the 
spin period is ~6 sec, and the angle between spin axis 
and sun-satellite line is always in the range 90° ± 20° . 
The positive ion detector (consisting of a hemispherical 
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electrostatic deflector and a Faraday cup) points perpen- 
dicular to the spin axis. Each measurement consists 
essentially of the positive ion flux integrated over an 
angle of 180° centered about the sun-satellite direction. 
Flux measurements are performed in 28 adjacent energy 
channels according to a programmed sequence, covering 
the range of energy per unit charge from about 200 to 
about 16,000 eV/Z. A complete energy distribution is 
obtained every 6.4 min, which is the duration of a 
telemetry subframe. The particular time sequence 
chosen for the energy channels allows a discrete-step 
coverage of the whole energy range in four subcycles of 
1.6-min duration. The sensitivity of the instrument 
is approximately 10 7 protons/cm 2 s or 5X10 6 a-par- 
ticles/cm 2 s. 

The magnetic field experiment on board HEOS 1 has 
been described by Hedgecock [1970] . 

AVERAGE PROPERTIES 

During the time period considered, a-particles were 
detected in 4,997 subframes, corresponding to a total 
observation time of ~20.8 days. The method of analysis 
to obtain the parameters of interest (a-particle number 
density N a> bulk velocity V a , most probable thermal 
speed Wy) has been described by Bonetti et al [1969] 
and Formisano et al [19706] . 
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Figure 1 shows the frequency distribution of individ- 
ual values of VjVp (ratio of a-particle to proton bulk 
velocity). The average value of VjVp is 1.02, identical 
to the value obtained by Robbins et al [1970] from 
Vela 3 observations. Though the deviation of this value 



Figure 1. Frequency distribution of individual values 
ofVjVp. 

from 1 is relatively small, attention should be paid to the 
lack of symmetry of the histogram. In several cases is 
definitely larger than Vp (in more than 300 spectra the 
difference is greater than 10 percent), while only a few 
spectra give V^jVp <0.9. It is interesting to note that 
the frequency distribution of values of VJVp given by 
Robbins et al [1970] exhibits the same asymmetry as 
the HEOS 1 histogram. On the other hand, theories of 
the solar corona expansion [Geiss et al , 1970; Nakada, 
1970] do not suggest any mechanism of preferential 
acceleration of helium relative to the hydrogen com- 
ponent of the solar wind (see also the discussion by 
Geiss, p. 566). 

The only experimental bias that could affect both the 
Vela and the HEOS histogram seems to be that when 
V a < V pi detection of a-particles could sometimes be 
prevented by the higher proton fluxes; indeed, if 
V a ~ 0.7 Vp, the a-particle peak in the energy per unit 
charge spectrum becomes coincident with the proton 


maximum and consequently cannot be observed. To 
evaluate the importance of the bias, we computed 
for each individual energy spectrum the lowest a-particle 
bulk velocity V am i n that still could be resolved from 
the proton spectrum (that is, which would give an 
observable secondary peak) when no second peak for 
a-particles was observed. Such computations were 
performed assuming W a = Wp and different values for 
the a-particle abundance (N a /Np). Results are shown in 
figure 2 in the form of frequency distributions of values 




Figure 2. Frequency distributions of values of 
Varmint Vp f or two a ~particle abundances . V a , m i n is 
defined as the lowest a-particle bulk velocity which still 
could be resolved from protons (i.e., which would give 
an observable secondary peak in the energy per unit 
charge spectrum). 

of V a m i n /Vp for two values of the ratio N^Np. As 
expected the bias decreases for high helium abundance, 
but it is certainly relevant for typical values of the ratio 
NJNp (~0.05). The behavior of the average values of 
V a and Vp versus N a (fig. 3) confirms this result; in fact 
for N a = 0.4 cm” 3 , V a reaches Vp , becoming smaller 
than Vp for higher helium densities. It should be jioted, 
however, that only a few spectra give V a < Vp. In 
figure 3 only cases with W a — Fp l > 20 km/sec are con- 
sidered consistent with detectability of a measurable 
difference between the two bulk velocities. 

We can conclude that within the limitation of present 
experiments, no difference between the average bulk 
velocities of the two species is observed; however, an 
average V a lower than Vp would have not been 
observed. Deviations of the ratio V^Vp from 1 appear 
to be present during particular time periods. 
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Figure 3. Average values of V a and Vp for N a larger 
than given values . Only cases with \V a -Vp\>20 kmjsec 
have been considered . 


The frequency distribution of the relative abundances 
of helium is given in figure 4. The average value of the 
ratio NoJNp is 0.051. Figure 5 shows a time history of 
the helium abundance with the measured sunspot 
number taken as an index of solar activity. Values of 
NJNp obtained by Vela 3 [Robbins et al , 1970] and 
by Explorer 34 [Ogilvie and Wilkerson, 1969] are shown 
with HEOS 1 averages for comparison. Averages from 
HEOS 1 and Vela 3 satellites have been performed over 



Figure 4. Frequency distribution of individual values 
ofNJNp. 


Figure 5. Time history of the ratio NJNp, according 
to observations performed by the satellites Vela 3, 
Explorer 34 and HEOS-1. A verages from HEOS-1 and 
Vela 3 [Robbins et al, 1970] have been performed over 
solar rotations (points) and over longer periods of time 
(bars). For Explorer 34 [Ogilvie and Wilkerson , 1969] 
only a total average is shown; the two values correspond 
to the average uncorrected for the instrumental bias 
(upper bar) and to the corrected average (lower bar). 
The number of cases over which the averages were 
performed is shown in the upper part of the figure. At 
the bottom , the measured sunspot number is given 
[Solar Geophysical Data , 1970]. 


solar rotations (points) and over longer periods of time 
(bars). For Explorer 34 only a total average was avail- 
able; the two values given in the figure correspond to the 
average uncorrected value (lower bar). The general 
pattern of the ratio N^/Np appears to confirm a 
dependence of helium abundance on the solar cycle. 

The plasma bulk velocity appears correlated with the 
helium abundance. Figure 6 shows the average values of 
F^and V a computed over the spectra with N^Np larger 
than 0.00, 0.01 , 0.02, .... For the study of the steady- 
state condition of the solar wind, data referring to 
postshock periods (defined as 24 hr following a sudden 
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Figure 6. Average values of Vp and V a computed over 
the spectra with N^Np larger than given values. The 
number of cases , over which the averages were per- 
formed, is shown on the top . 

commencement) have been removed from these aver- 
ages. An increase of ~120 km /see is observed in the bulk 
velocity of both species when NJNp increases. 

Figure 7 presents a comparison between the Jjehavior 
of V a , Vp versus N a and Vp versus Np. The V p versus 
Np plot, which was made with all the proton data 
(fig. 7(b)), shows the well-known result that the proton 
flux NpVp is approximately constant, in contrast to the 
general pattern of the V a versus N a plot (fig. 7(a)), 
which shows a tendency to an increase of V a with N a . 
Again it should be noted that for high densities of 
a-particles, V a becomes lower than Vp. 

number of cases number of cases 




Figure 7. Average values of V a and Vp plotted versus 
N a (a) and of Vp versus Np (b). Number of cases is given 
on top. The Vp-Np plot was made using all the available 
data. 


The increase in helium abundance together with the 
increase in the bulk velocity is predicted by the diffusion 
models of the solar corona [Nakada, 1970]. An increase 
in the coronal temperature, indeed, should increase the 
solar wind bulk velocity following the magnetohydro- 
dynamic models [ Parker , 1963; Hartle and S tarrock, 
1968] and should also decrease the diffusion time, 
increasing the helium abundance in the expanding 
corona. 

TWO-BEAM INSTABILITY 

As shown in the previous section and by Formisano et 
al. [ 1 970Z?] it is possible, occasionally, to observe 
different bulk velocities for protons and a-particles. 
Theoretically, it has been predicted [Geiss et al., 1970; 
Nakada , 1970; Yeh, 1970; AUoucherie, 1970] that the 
helium bulk velocity should be lower (about 20 percent 
or more) than the proton bulk velocity in the solar 
corona, while some instability of the two stream type 
should be operating in the solar wind to equalize Vp and 
V a - The problem has been studied with our data. The 
observations should reveal the quasi-steady state quickly 
reached by the solar wind at a few tenths of AU from 
the sun. Figure 8 shows the behavior of Wp, W a versus 
VJVp. The protons and helium average most probable 
thermal speeds increase for values of V a and V p 
different from 1 . The observed variations of W a and Wp 
correspond to an increase of the proton and a-particle 
temperature of a factor 2. This result was already shown 
by Formisano et al. [19706] . Figure 8 refers to more 
data, and the result is confirmed. 



Figure 8. A verage values of W a and W p plotted versus 
VoJVp. Averages have been performed on VjVp inter- 
vals of 0.05. In the upper part of the figure histograms 
of the total number of cases considered in the averages 
are given. 
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The two-stream instability predicts that the faster 
beam will be slowed down, giving energy to waves and 
increasing the thermal energy of the particles. If the 
faster beam is the major beam— that is, has higher 
density (protons in the solar wind)— the energy available 
for waves and thermalization is large, and therefore a 
more rapid increase of the temperature of the slower 
beam and a higher noise level of unstable waves is 
expected. If the faster beam is the minor one (lower 
density, a-particles in the solar wind), a slower increase 
of the thermal energy of the slower J>eam is predicted. 
This behavior is shown in figure 8; W a increases faster 
(for V a < V p ) than W p (for V a > V p ). 

Figure 9 shows the actual solar wind points in relation 
to the two stream instability criterium of high frequency 
longitudinal waves for two maxwellian distributed beams 
[ Taylor ; 1970], In a plane b y c [b - 4 (WqJW p ) 2 , 
c= W a -V p \ fW p ) plasma will be stable or unstable 
depending on the quantity a = (N^Np) (W p /W a ) 3 . The 
stable region is the one on the left of the line corre- 
sponding to the given value of the parameter a. 

Figure 9 shows for each (b 9 c) pair the edge of the 
stable region. The “a” scale is given at the top of the 
figure. All the data plotted have W a ^V p \> 10 km/sec. 
All the points fall within the stable region, some just 



Figure 9. The scatter plot of b = (Wp/W#) 2 and 
c = Wq- V p \/W p . The two-stream instability (for longi- 
tudinal waves) lines are shown for different values of the 
parameter a ~ Ny/Np (W p /W a ). Black areas represent 
the edges of the stable region for each (b,c) pair. 


behind the stability line, most well inside the stable 
region. Lower values of \V a - V p \ would give more 
points in the more highly stable region. 

It should be remembered, however, that the two- 
streaming plasma is unstable with respect to both 
longitudinal and transverse waves. Parker [1961] has 
shown that Alfven waves could be generated by a 
two -streaming plasma. Kennel and Petschek [1968] 
have related the firehose instability to the two-beam 
instability and the resulting instability criterium 
becomes 


where 



P,e,a 


3 mW 2 

4 B 2 /8tt 


( 1 ) 


is the ratio between the total particle pressure and the 
magnetic pressure, Oi p (oi a ) is the ratio between the 
proton (a-particle) number density and the total ion 
number density; A/3 = j3 1 | - j3 1 | = 3/3 - 2j3j_; and 
/3|l/j3 i~K is the anisotropy ratio. It is evident from 
equation (1) that the firehose instability is modified if 
a-particles have V a * Vp and that the difference 
V a —Vp must be reduced by a firehose stable anisotropy 
in the three species. Figure 10 shows the stable region 
with respect to both firehose and two stream instability. 



Figure 10. Stability region with respect to the firehose 
and two stream instability . Equation (1) has been 
plotted for j 3a a a p ~ 0.11 and 0.58. The 0.11 value 
refers to the set of parameters given on the top of the 
figure (the corresponding position in the A/3, A V/W p 
plane is represented by the black point). 
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As an example, a possible set of parameters and the 
corresponding position in the Aj3, A V/Wp plane have 
been indicated. 

It should also be noted that another instability, of the 
resonant type, has been found by Barnes , [1970] for 
bow shock reflected protons. It is not yet clear whether 
this instability can be important for the a-particle and 
proton streams. 


SHOCK WAVE OBSERVATIONS 

The behavior of a-particles through two interplanetary 
shock waves (February 28 and March 19, 1969) has been 
studied by Formisano et al [1970/?] using the HEOS 1 
plasma data. Here we present the a-particle observations, 
together with interplanetary magnetic field measure- 
ments, during the two above-mentioned shocks and a 
third one, which occurred on March 25, 1969. This last 
event is discussed in detail by Chao et al [1971] . 

Figure 11 shows the event, which occurred on February 
28, 1969, detected at ground as a sudden commencement 
at UT 0425, by 19 stations. At UT 0422 a sudden increase 
of a-particle density, bulk velocity, and thermal speed is 
observed; V a changes from 507 km/sec to 637 km/sec; 
N a from 0.35 cm” 3 to 0.80 cm” 3 ; W a from less than 
57 km /sec to ^88 km /sec. 

The proton behavior is quite different: during about 
half an hour (from UT 0405 to UT 0431) Vp changes 
only ~50 km/sec. The proton density, which has an 
increment of ~30 percent at UT 0412, decreases to the 
previous value of 5 cm” 3 when the a-particle disconti- 
nuity is detected. A small gradual increase is observed 
for the proton thermal speed. The proton energy 
spectrum, which was previously maxwellian, exhibits a 
non-maxwellian high-energy tail after the a-particle 
discontinuity, as had been shown by Formisano et al 
[19706] . The a-particle parameters show the character- 
istic variations of a shock. The identification of the 
a-particle discontinuity with a shock wave is confirmed 
by the magnetic field measurements. Indeed, simultane- 
ously with the a-particle discontinuity, the magnetic 
field intensity increases of a factor 2 (from —67 to 
~ 12 7 ). 

Magnetic field oscillations are noted in the figure, 
before the sudden variation at UT 0422. They appear to 
be related to the earth’s bow shock crossed by the 
satellite, moving outward from the earth, at UT 0120; 
during the 3 -hr period from UT 0120 to UT 0422, these 
oscillations are permanently detected by HEOS 1 , indi- 
cating the presence of waves associated with protons 
reflected from the earth’s bow shock [Asbridge et al , 
1968]. 



Figure 1 1 . Positive ion parameters and magnetic field 
for the event of February 28, 1969. From the top are 
given: the magnetic field intensity , the a-particle to 
proton kinetic energy density ratio; the a-particle and 
proton number densities; the most probable thermal 
speeds of protons and a -particles; the proton and 
a-particle bulk velocities. Broken lines refer to 
a-particles, continuous lines to protons. 

Figure 12 illustrates the March 19 event observed at 
ground as a sudden commencement at UT 1959. Three 
minutes before (at UT 1956 with an error of ±0.8 min) 
the protons exhibit a velocity discontinuity of 
~45 km/sec (from 396 km/sec to 439 km /sec) and a 
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Figure 12. Positive ion parameters and magnetic field 
intensity for the event of March 19, 1969 . 

density, increase by a factor of 4 (from 7 protons/ cm 3 to 
27 protons/cm 3 ), while the thermal speed increases by 
only 10 km/sec. The a-particle discontinuity ( V a from 
404 km /sec to 453 km /sec and N a from 0.2 cm" 3 to 
0.5 cm" 3 ) is observed with a time delay of —3 min 


compared to the proton discontinuity [Formisano et al, 
1970&). The magnetic field measurements show a 
sudden change of the field intensity (from ~5y to ~107) 
at 1956 UT, simultaneously to the proton discontinuity. 
Later on, the magnetic field undergoes strong fluctua- 
tions. 

Figure 13 illustrates the third shock wave occurring on 



DAY 84,1969 

Figure 13. Positive ion parameters and magnetic field 
intensity for the event of March 25, 1969 . In the 
number density plot the scale of N a should be multi- 
plied by KT 1 . 


594 



March 25. At UT 0646 a sudden variation of the proton 
parameters is observed ( V p from 435 km/sec to 
490 km/sec; N p from ~6 cm“^ to 10.4 cm” 3 ; W p from 
40 km/sec to 60 km/sec); at the same time, the magnetic 
field intensity increases from IO7 to I87. The a particles 
exhibit a different behavior. No relevant change of V a is 
observed when the shock is detected, while a large 
discontinuity is observed a few minutes before the 
shock: V a goes from 430 km /sec at UT0614 to 
515 km/sec at UT 0633 (two measurements are missed 
between). 

Due to the data gap, we have no information on N a at 
UT 0633, when the discontinuity of the a-particle bulk 
velocity is observed; however, the fact that a particles 
escape detection in the two subframes between UT 0614 
and UT 0633 may suggest that N a was lower during that 
period. When the proton shock is detected, at UT 0646, 
N a increases from 0.3 cm -3 to 0.7 cm” 3 . 

The behavior of protons and a particles through the 
shock is better illustrated in figure 14. Three energy per 
unit charge spectra are shown: subframe 25630 is 
observed before the V a change; subframe 25633 is 
observed just after the V a increase, but before the 
proton discontinuity; subframe 25635 is taken just 
behind the shock. A maxwellian distribution was fitted 
to the three highest fluxes for both protons and 
a-particles; the obtained parameters are shown for both 
species. In all cases the maxwellian fits do not appear to 
be very good; in particular, a non-maxwellian high 
energy tail is clearly observed at UT 0633, when there 
was a large difference between proton and a-particle 
bulk velocities. 


CONCLUSIONS 

We can summarize our results as follows: 

1. Within the limitations of present experiments, no 
difference between the average bulk velocities of 
protons and a-particles is detected. However, due to 
an instrumental bias, the possibility of V a being lower 
than Vp cannot be ruled out. 

2. Data from the satellites Vela 3, Explorer 34, and 
HEOS 1 give evidence of a dependence of helium 
abundance on the solar cycle. 

3. On average, high helium abundances are observed 
together with high bulk velocities, in agreement with 
diffusion models of the solar corona. 

4. When a definite difference between proton and 
a-particle bulk velocity is observed, the stability with 
respect to high frequency longitudinal waves is 
ensured by higher thermal speeds. However, the 
problem of transverse wave instability, due to differ- 
ences between Va and V p , has not been investigated. 


MARCH 25.1969 

• experimental points 
x fit with given parameters 



Figure 14. Positive ion energy spectra observed on 
March 25 , 1969 , at UT 0614, 0633 and 0646. Electrom- 
eter counts (black points) are plotted versus energy per 
unit charge . P indicates proton peaks; a indicates 
a-particle peaks. Crosses indicate the maxwellian fit of 
the observed spectra for the given values of proton and 
a-particle parameters . 

5. The behavior of a-particles through interplanetary 
shock waves is not a simple one. Strong disconti- 
nuities may be observed in one species, without 
occurring in the other, or having different character- 
istics. In particular, the magnetic field observations 
confirm the interpretation of the a-particle disconti- 
nuity, occurred on February 28, 1969, as a shock 
wave. 

6. When a strong discontinuity is detected only for 
a-particles, the proton distribution function shows a 
large non-maxwellian high-energy tail. 
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DISCUSSION J. Hirshberg Remarking on the helium abundance increasing as the solar cycle goes 
on, I compared your data with the Vela data and it seems that in those two samples the 
most probable value remained stationary. What happened was that as the solar cycle 
progressed you obtained more high helium observations which would be consistent with 
the notion that enhanced helium was coming out of flares. 

G. Moreno That’s right, the modal value was the same in two cases, only the 
enhanced helium changed. I forgot to say that the correlation of the bulk velocity of both 
protons and a particles bulk velocity with the helium abundance has been done only 
for quiet periods. 

M. Dryer How many energy channels do you have in your analyzer? 

G. Moreno Twenty-eight. Computation suggests that the velocity spread from one 
channel to the next error in velocity is perhaps less than 20 km; in many cases when 
a particles are in three channels it becomes much less. 

D. Heymann Several years ago, as you well know, it was suggested by Dr. Michel that a 
tenuous atmosphere on the moon is removed by the interaction of this atmosphere with 
the solar wind rather than by gravitational escape. 40 A was not removed from the 
moon but was punched back into the lunar surface by essentially the same mechanism 
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that Michel proposed for the escape of neutral atoms, i.e., via ionization then acceleration 
in the interplanetary electric and magnetic field from the moon. The 3 6 A, as he argued, 
varies quite a bit, but it varies because it is surface correlated. He argued and showed 
that all the 36 A is essentially from the solar wind, and implanted in the surface 
of these particles. Now, the curious thing is that when we measured 40 A in the same 
particles they correlate, I mean, they co-vary with 3 6 A. And in fact 40 A/ 3 6 A in these 
samples is about unity; that is a paradox because all theoreticians tell us that 40 A/ 3 6 A 
in the sun ought to be much less than unity, in fact perhaps much lower than 10~ 4 . So 
the question is where the 40 A that we are seeing in the samples comes from. It cannot 
come from the solar wind proper. You may wonder and say perhaps this 40 A has arisen 
from 40 K decay in these samples and what we note is simply that these samples are low 
in 40 K or have a young age or both, and that these samples have much potassium or have 
an old age or both. But that cannot be true because for any of these points if you 
calculate an honest 40 K/ 40 A age you find that all the ages are greater than the accepted 
age of the moon and in fact calculate ages greater than 7X10 9 yr. So the 40 A Has been 
produced from potassium not in these samples but somewhere else, has been separated, 
and then reimplanted into the lunar regolith. In other words, what we are looking at is 
Argon which was produced in the body of the moon, by potassium decay, was made 
available in the lunar atmosphere. You can, for example, envisage Argon just oozing out 
slowly. It might be that there were large impacts on the surface of the moon which threw 
out vast amounts of rock and that Argon was subsequently reimplanted— being coupled to 
the solar wind which is the driving agent of this reimplantation into the soil— and that is 
basically the reason we see the covariation of 4 0 A and 3 6 A. 



